Kinetics of the evaporative cooling of an atomic beam 
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We compare two distincts models of evaporative cooling of a magnetically guided atomic beam: 
a continuous one, consisting in approximating the atomic distribution function by a truncated 
equilibrium distribution, and a discrete-step one, in which the evaporation process is described in 
terms of successive steps consisting in a truncation of the distribution followed by rethermalization. 
Calculations are performed for the semi-linear potential relevant for experiments. We show that it 
is possible to map one model onto the other, allowing us to infer, for the discrete-step model, the 
rethermalization kinetics, which turns out to be strongly dependent upon the shape of the confining 
potential. 

PACS numbers: 32.80.Pj,03.75.Pp 



I. INTRODUCTION 

Evaporative cooling is a very powerful technique 
that allowed the achievement of quantum degeneracy in 
dilute atomic vapors 2|. On the theoretical side, apart 
from direct numerical simulations Q, several models of 
evaporative cooling have been studied, which can be clas- 
sified in two categories: continuous and discrete ones. 

In discrete models, the process of evaporative cooling is 
approximated as a series of truncations of the atomic dis- 
tribution function, followed by rethermalization towards 
equilibrium 4]. Through many steps, the phase-space 
density of the atomic sample increases up to order unity. 
The advantage of such models lies in their simplicity. 
However, in the case of trapped clouds of atoms, those 
models are not realistic, since, experimentally, the evap- 
oration is done by ramping down continuously a radio- 
frequency knife. Moreover, such models give no indica- 
tion on the kinetics of the evaporation. Therefore, one 
needs to resort to more elaborated, continuous models, 
using an appropriate Ansatz for the distribution function, 
namely a truncated equilibrium function. It is then pos- 
sible, starting from Boltzmann equation, to obtain the 
time evolution of the temperature and of the number of 
atoms . Quantitative comparisons between the predic- 
tions of those two types of models have not been made 
to date. 

Evaporative cooling has been revisited in the context 
of the cooling of a guided atomic beam, in view of achiev- 
ing a continuous-wave atom laser. The proposal |(j used 
a continuous evaporation model to predict the possibil- 
ity of achieving quantum degeneracy by using transverse 
evaporation on an atomic beam, confined transversally 
by an harmonic potential. This proposal triggered ex- 
perimental work in several groups, and a discrete model 
of the evaporative cooling of a beam was developed in 
close connection with recent experiments 0, || . 

In this paper, we address the problem of the trans- 
verse evaporative cooling of a magnetically guided atomic 
beam with an energy-selective "knife". We first de- 
velop a continuous model using hydrodynamic equations, 
adapted from 6], for an experimentally realistic trans- 



verse potential. Then we describe the same process with 
a discrete-step evaporation model, analog to the one used 
in [3 • Finally, we compare the results given by those two 
distincts models. In particular, this comparison allows 
us to study the influence of the shape of the confining 
potential on the kinetics of rethermalization and on the 
evaporation 'ramp'. 

We consider an atomic beam with a flux <I>, a mean 
longitudinal velocity v, and a temperature T, propagat- 
ing in a quadrupole magnetic guide of axis z, providing 
the following semi-linear transverse potential : 



U(r) = fi^B 2 + b 2 r 2 _ M B Q 



(1) 



Here, fi is the magnetic moment of the atoms, b the 
transverse gradient of the two-dimensional quadrupole 
magnetic field created by the guide, r = x 2 + y 2 the 
distance from the guide axis, and Bo a longitudinal bias 
field used to avoid Majorana spin-flips @- The on-axis 
potential is taken as the origin of energies. One defines 
the dimensionless parameter a = fiBo/k^T. For typi- 
cal experimental parameters, the evaporation starts with 
a <C 1 where the potential experienced by the atoms 
is essentially linear, and degeneracy is reached in the 
regime a 3> 1, where the potential is essentially har- 
monic. Therefore it is crucial to take into account the 
real shape of the potential in order to describe the whole 
evaporation process. 

An important quantity characterizing the beam is the 
s-wave elastic collision rate 7, given by: 
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(2) 



where a is the s-wave scattering cross-section. 

Two-dimensional transverse evaporation is applied 
in order to selectively remove atoms from the beam. 
In practice, this is achieved by driving transitions to 
an untrapped state with radio-frequency or microwave 
fields 0. The evaporation criterion then relates to the 
transverse energy e and to the angular momentum of the 
atom along z. For the sake of simplicity, as in Ref. 
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we assume in the following that any atom having a trans- 
verse energy e r/k^T, where 77 is the evaporation pa- 
rameter, is evaporated. This criterion implicitly assumes 
sufficient ergodicity of the atomic motion. Experimen- 
tally, this simple energy criterion can be well approxi- 
mated by the multi-radii evaporation scheme described 
in0. 



II. CONTINUOUS MODEL 

In this section, we assume that the evaporation takes 
place over the whole guide length. The height of the 
energy knife £ e v(z) = f]kg,T depends on z in order to 
perform forced evaporation. We therefore assume that 
the beam's distribution function is a local equilibrium 
one p|, truncated at the energy e GV . 

By using such an Ansatz in the Boltzmann equation, 
one gets a set of coupled hydrodynamic equations relating 
the following local quantities characterizing the beam: 
the linear density n(z), the longitudinal velocity v{z), 
and the temperature T(z). The details of this somehow 
lengthy calculation can be found in ref. 0, for the 
case of a harmonic transverse potential. The validity of 
such an Ansatz was checked with a molecular dynamics 
simulation of the process. In the following, we present 
the generalization of this analytical approach to the case 
of a semi- linear confinement. 

In the stationary regime, the hydrodynamic equations 
read: 




Temperature T (jiK) 

FIG. 1: Evaporation "trajectories" in the temperature/flux 
plane, for evaporation with a constant 77 for the continuous 
model (solid lines) and for the discrete one (dots, each symbol 
representing the effect of one evaporation step). The back- 
ground grey scale (with the white labels) shows the on-axis 
phase-space density. The qualitative shape of those trajec- 
tories are similar for both models. For n high enough (e.g., 
77 = 8) , the trajectories are very close in both models. The 
number of evaporation zones required for reaching degeneracy 
in the discrete-step model, for rj — 5 (resp. 6, 8), is 88 (resp. 
152, 526). 
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(3a) 
(3b) 
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where v t h — y feT/m. They correspond, respectively, to 
the evolution of the flux, of the momentum, and of the 
enthalpy of the beam. The notation (U) stands for the 
thermal average, at the temperature T, of the potential 
energy U. These equations are well suited to describe a 
supersonic beam with a high enough Mach number (typ- 
ically v > 3v t h)- 

In the semi-linear potential 0J, the mean value of the 
potential energy reads (U) = k^T(2 + a)/(l + a). T\ and 
T2 correspond to the evaporation-induced particle and 
energy loss rates, respectively. They are proportional to 
the elastic collision rate 7 and obviously depend on the 
evaporation parameter 77: 
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T i = 1 ^--K i (r,,a) (i = l,2). (4) 

The functions K i are given by the following integral: 
/■ 77+1/2 

Ki(f],a)= f(x,a)gi(x,T])dx, (5) 

Jo 



with f(x,a) = z 3 / 2 (5a + 2x)er r >- x l 2 j {\ + 2a) being a 
contribution from the density of states per unit length in 
the semi-linear potential, g\{x, rj) — e~ x (r] — x — 1/2) + 



e 1 1/2 , and g 2 (x,ri) 
l/2) 2 -2-3x/2-T?x]. 



-7,-1/2 



(3 + 2t7-.t) +e- x [{r ] + 



The solid lines in Fig. ^ depict the evolution of the 
beam's flux and temperature obtained by solving the hy- 
drodynamic equations assuming that the evaporation 
parameter 77 remains constant throughout the evapora- 
tion. The initial conditions are the experimental ones of 
Ref. 0, in which 87 Rb (p = /jb/2) atoms are used: at 
z = 0, one has $ = 7 x I0 9 at/s, v — 60 cm/s, and 
T = 570 fjK. The gradient is 800 G/cm and a B = 1 G 
bias field is applied. When one increases the value of 77, 
degeneracy (phase-space density p ~ 1) is achieved for 
higher fluxes (and therefore higher temperatures) since 
the evaporated particles are very energetic and conse- 
quently the evaporative cooling is more efficient. The 
change in the slope of the "evaporation trajectories" for 
T ~ nBo/kB ~ 30 /xK is due to the fact that the confine- 
ment experienced by the atoms changes from essentially 
linear (a = 0.06 <C 1) to essentially harmonic (a ^> I) as 
the temperature is reduced. Indeed, the gain in phase- 
space densityscales differently with the shape of the po- 
tential HH13- 
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III. DISCRETE-STEP MODEL 



We now turn to the description of the evaporative cool- 
ing process with a discrete-step model. One evaporation 
step reduces the atomic flux from $ to After rether- 
malization, the beam acquires a new temperature T' . In 
order to calculate the relative variations of flux <p = 
and of temperature r = T'/T, we adapt the approach of 
Ref. to the semi-linear potential and to the energy- 
dependant evaporation criterion. We therefore introduce 
the two-dimensional density of states q{e) in the semi- 
linear potential QJ: g(e) oc e[2 + e/(fiBo)]. The fraction 
if of remaining atoms after one evaporation step is: 



ip(r),a) 



J£ B g{e)e-c/ k * T &e 
r 3 00 e(e)e- £ / fcBT de 
_ 2 + 2a{\ + 77) + 77(2 + 77) ^ 
2(1 + a) 



(6) 
(7) 



In order to derive r, we first calculate the transverse en- 
ergy e of the remaining atoms: 



/ °° g(e) e- £ / kBT de 



(8) 



We define the dimensionless parameter #(77, a) = 
ej (/cbT) 01 • The conservation of the total energy during 
rethermalization gives: 



k B T ( (9$ 



,3fc B T' 



(9) 



where the average (U) T , , of the potential is taken at 
thermal equilibrium with a temperature T", i.e. with 
a' = cv/t. This yields a quadratic equation in r, with 
the solution: 



2(9-1- — 5ay5 + y/28pa(26 + tp) + (28 + ip- 5aip) 2 
T ~ lip ' 

One then readily obtains the relative variations of the 
collision rate 7 and of the phase-space density p after an 
evaporation step. 

The corresponding "evaporation trajectories" (for 77 
constant) are depicted with circles on Fig. ^ Each sym- 
bol represents the flux and temperatures ($ n , T n ) of the 
beam after the n th evaporation zone. 



IV. DISCUSSION 

We compare the results given by the two models, in 
terms of evaporation trajectories and of the efficiency of 
evaporation. As depicted on Fig. ^ the evaporation tra- 
jectories have the same qualitative behavior in both mod- 
els. However, for a given 77, the discrete-step evaporation 
leads to higher final fluxes and temperatures. For a high 



u 

a 

Efi 
W 



4 - 



FIG. 2: Evaporation efficiency £ as a function of the evap- 
oration parameter 77, for a harmonic transverse confinement. 
The solid (resp. dashed) line corresponds to the continuous 
(resp. discrete-step) evaporation model. 



evaporation parameter (77 = 8), the trajectories given by 
both models almost coincide. 

To make these statements more quantitative, we intro- 
duce the figure of merit for an evaporative cooling ramp, 
i.e. the relative variation of the beam's phase-space den- 
sity p for a given loss of particles. We therefore define 
the evaporation efficiency £ as: 



A In p 
dln<f> 



(10) 



This quantity is straightforward to calculate for a given 
evaporation model. Fig. represents £(77) for the case 
of a harmonic transverse confinement. As expected, £ 
increases with 77 in both models. It appears that discrete- 
step evaporation is more efficient than the continuous 
one, which can be understood qualitatively by the fact 
that in the latter scheme, some atoms are evaporated 
without giving rise to a temperature reduction, a process 
commonly called "spilling" PIUI- However, for r\ high 
enough, the efficiencies of both models almost coincide. 

We now turn to the kinetics of evaporation: for a 
given shape of the potential (linear or harmonic) and 
a given value 77D of the evaporation parameter in the 
discrete model, we determine the corresponding parame- 
ter t\c(j]d) for the continuous model, which leads to the 
same 'evaporation trajectory' in the (T, $) plane. This 
mapping allows us to extract information on the kinet- 
ics aspects of the rethermalization between evaporation 
zones for the discrete models. More precisely, we infer the 
number N c of collisions required to rethermalize between 
successive zones. For this purpose, we integrate over 
time the collision rate from the continuous model, with 
an evaporation parameter r/ci^D), between two points 
(T n ,$ n ) and (T n+ i, $ n +i) of the evaporation trajectory 
obtained with the discrete model. 

To allow for a quantitative comparison between lin- 
ear and harmonic confinements, we scale the evapora- 
tion parameter 770 as in Ref. Q by defining iju — i]d/2 
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FIG. 3: Number N c of collisions necessary for rethermaliza- 
tion between two evaporation zones (inferred by mapping the 
continuous evaporation model onto the discrete one, see text) 
as a function of the scaled evaporation parameter f\o- The 
solid (resp. open) circles correspond to a linear (resp. har- 
monic) confinement. 



(resp. f\u = ?7_d/3) for a harmonic (resp. linear) confine- 
ment. Evaporation with a given normalized parameter 
t)d then yields approximately the same flux reduction, 
independently of the shape of the confining potential. 
On Fig. we have depicted N c as a function of t)d , for 
both linear and harmonic confinements. As expected, N c 
decreases with increasing f)D, since the atomic distribu- 
tion is less and less affected by the evaporation, yield- 
ing a faster relaxation towards equilibrium. For a linear 
confinement, 50% to 100% more collisions (as compared 
to the harmonic one) are required for rethermalization 
between evaporation zones. This dependence of the ki- 
netics on the shape of the confining potential is remi- 
niscent of what is known for thermalization of confined 
gas mixtures 0] : the rethermalization time is shorter in 
a homogeneous system than for a trapped cloud. In a 
power-law trap of exponent 5, the rethermalization time 
decreases when S increases, which simply originates from 
the different scaling laws of the density of states. 

Therefore, two competing effects need to be considered 
when one studies the whole evaporation process: in a lin- 
ear potential, the kinetics is slow but the gains in collision 
rate and in phase-space density scale more favorably Q 
than in a harmonic confinement. In terms of the min- 
imum number of collisions required to achieve a given 
gain in phase-space density, those two effects compensate 
each other. For instance, we find that at least 500 (resp. 
630) collisions are necessary to gain a factor 5 x 10 7 in 
phase-space density, in a purely harmonic (resp. linear) 
potential, for an evaporation parameter rjr> ~ 4.5 (resp. 
r\Y) ~ 5.5). However, in terms of evaporation length, the 
difference between harmonic and linear confinements is 
still large, as runaway evaporation can only occur in the 
latter case for a two-dimensional evaporation. As an ex- 
ample, we consider two beams with the same initial flux 
(7 x 10 9 atoms per second) and temperature (200 /jK), 




Initial flux $i (10 9 at/s) 

FIG. 4: Minimal evaporation length L ov needed to reach a 
phase-space density p = 1 with evaporation at constant n for 
an atomic beam of initial flux $i and temperature Tj, propa- 
gating at 60 cm/s in a guide with a gradient b — 800 G/cm 
and a longitudinal bias field Bq = 1 G. The iso-contours of 
Lev show the large advantage of starting with a temperature 
on the order of 100 /iK for reaching degeneracy in a reason- 
able length. The optimal evaporation parameter r\ is almost 
constant (~ 6) over the whole parameter space explored here. 

with an initial phase-space density p\ ~ 8 x 10~ 7 , prop- 
agating at 60 cm/s either in a purely harmonic guide, 
either in a purely linear one. For the former case, the 
initial collision rate is 7 ~ 37 s _1 , and slightly decreases 
to 7 ~ 29 s -1 after evaporation to degeneracy. The evap- 
oration is performed at rjc = 6, a value that minimizes 
the evaporation length L cv , which reaches about 11 m. 
For a linear confinement, although the initial collision 
rate is only 7 ~ 19 s" 1 , its final value reaches 380 s _1 
due to the runaway character of the evaporation. The 
total length needed is only 6 m. Interestingly, the total 
number of collisions that actually occurred within the 
beam is almost the same for both confinements (~ 550), 
as is the evaporation parameter minimizing L ev . 

In practice, due to Majorana spin-flips |9j, the guide 
potential needs to be semi-linear. For a guide with pa- 
rameters b = 800 G/cm and Bq = 1 G, we have studied 
the minimal evaporation length L cv needed to reach de- 
generacy, as a function of the initial flux and initial 
temperature 2j, assuming a beam velocity of 60 cm/s. 
The result is plotted on Fig. and shows as expected 
that Lev decreases with lower 1] and higher <I>j. The 
optimal evaporation parameter is almost constant with 
the value 77 = 6 for our range of parameters (<I>i, T|). The 
evaporation length determined this way is very well fitted 
by a function of the form 

j,3/2 

L ov =S Lo-^- • (11) 

This scaling can be easily understood: since here the 
runaway effect exists only at the very beginning of the 
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evaporation ramp, before the effective shape of the poten- 
tial becomes harmonic, the evaporation length is simply 
inversely proportional to the initial collision rate. For 
an initial flux 7 x 10 9 atoms/s, an initial temperature 
200 /iK, a guide gradient b = 800 G/cm and a bias field 
Bq = 1 G, quantum degeneracy is reached for L cv ~ 7 m, 
which shows the beneficial influence of the increase of the 
collision rate at the beginning of the evaporation ramp 
when the potential is essentially linear. 

Actually, the evaporation length deduced here could be 
reduced by lowering the beam's mean velocity (e.g, with 
the use of a tilted guide as it cools down, provided 



that the beam stays supersonic enough. 
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